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Abstract 

A large class of stellar systems (e.g., planetary nebulae 
(PNe), supernova envelopes, LBV stars, young stars in 
formation) shows structures in their accretion/ejection 
phase that have similar characteristics. In particular, 
one currently observes for these objects equatorial 
discs, axial ejections and stable bipolar shells. How- 
ever these simple shapes, which are expected to be 
solutions of standard hydrodynamical equations, are 
not yet fully understood. In this paper, we suggest a 
new form of these equations that takes into account 
the fractality and the irreversibility of particle motion 
in such processes. Then we study in this framework a 
general infall or ejection motion in a central spherical 
potential. From the stationary solutions allowed by 
this new hydrodynamical system, we deduce a specific 
distribution of matter density, described in terms 
of probability density for ejection angles. A global 
classification of predicted shapes, depending on the 
values of conservative quantities such as (E 2 ,L 2 ,L Z ), 
is given. These results are compared with the available 
observational data, and allows us to theoretically 
predict the possible existence of more complicated 
structures and of correlations between observable vari- 
ables, which could be checked by future observations. 

keywords: relativity - gravitation - hydrody- 
namic - planetary nebulae - LBV star - supernovae - 
star formation 



1 Introduction 

With the recent evolution of observational techniques, 
a great diversity of unexpected gravitational structures 
has spring up in many astrophysical domains. The 
large number and the stability of these structures has 
pointed out the incompleteness of their understanding 
in terms of standard approaches. This still represents 
a large and fundamental open problem. Moreover, one 
observes similar structures in many different systems, 



at many different scales and in various conditions char- 
acterizing the underlying medium, such as stable discs, 
rings, disc/jets combinations and bipolar outflows. 
The present works on the formation and the dynamics 
of planetary nebulae are based on hydrodynamical col- 
lisions between slow and fast winds emitted at different 
stellar evolution stages. This is the Interacting Stel- 
lar Wind (ISW) model, generalized by Balick 0. In 
this model, slow winds (10-20 km/s) are firstly emitted 
during the AGB phase then, at the end of this period, 
UV-radiations ionize this circumstellar slow medium. 
After this initial ejection, the stellar surface still losses 
its matter, but fast winds (sa 1000 km/s) are observed 
and enter in collision with slow winds. This collision 
induces a compressed zone where a hot bubble creates 
the observable bright shells of planetary nebulae. The 
spherical and elliptic shapes are satisfactorily under- 
stood in the framework of this description ([Jj], |lfi|L 
but the majority of the planetary nebula shapes cannot 
result from this model |13| . Despite the improvement 
of hydrodynamic simulations, even the simplest bipo- 
lar shells are not found to be formed in an universal 
way , since each particular case needs an ad-hoc hy- 
pothesis on the density of presumed slow equatorial 
winds and observed fast winds. Moreover, bright bipo- 
lar shells are not directly connected to the hot bubble 
resulting from the collision Thus, in the standard 
hydrodynamical approach, the initial hypothesis made 
about the nature of the slow and fast winds (which are 
connected to each central system studied) play a cen- 
tral role. 

In the present paper, we shall explore a different though 
complementary approach to these problems. Indepen- 
dently of the initial wind densities, we can study the 
consequences of the specificity of individual trajecto- 
ries of particles making the shells. We use a fluid- 
like description, in which the possible trajectories are 
non-deterministic. Because of the information loss due 
to collisions and to the continual ionization/diffusion 
of central star photons, individual particle trajectories 
can no longer be strictly followed, so that a statistical 
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description in terms of fractal and non-differentiablc 
trajectories is adopted. Such a description finally gives 
a new form to the hydrodynamical equations. 

2 Theory 

In the following theoretical development, we use gen- 
eral hydrodynamical concepts, then we apply them to 
the specific case of an ejection process. It can be al- 
ready noted that these solutions can be also used, with 
a mere time reversal, for describing infall motion in a 
spherical potential. 

In a typical process of stellar mass ejection (developed 
above), the interaction between winds or ISM and the 
continual interaction of emitted photons allow one to 
set the following simple hypotheses: 

(i) For each particle, the information about its tra- 
jectory is lost beyond some small space-scale and time- 
scale (compared with the global space-scale and time- 
scale of structuration of the system). Such a memory 
loss of the previous dynamical conditions is compati- 
ble with a fractal description using a fractal dimension 
Dp = 2 for trajectories (as in standard Brownian mo- 
tion). Only global informations are conserved. 

(ii) The test particles can follow an infinity of poten- 
tial trajectories, and the study of such systems there- 
fore comes under a fluid approach (like in their stan- 
dard description). 

(iii) The last main point is the irreversibility of 
this process. This irreversibility is local and global: 
for small scales, this is connected with the non- 
deterministic trajectories, and at the scale of the global 
system, this represents the impossibility for the system 
to come back to its initial conditions. 

These three points can be summed up in terms of three 
fundamental conditions which lead to the construction 
of new hydrodynamic tools: 

(i) The fractality of each individual trajectory be- 
yond some transition scale. It is described in terms of 
elementary displacements which read: 



dX = dx + d£, 



(1) 



where dx = vdt is a classical, differentiable vari- 
able, while dt; is a non-differentiable, stochastic vari- 
able, that describes the fractal fluctuation. Its non- 
deterministic character implies that it is known only 
through a statistical description. Namely, in the sim- 
plest case only considered here (fractal dimension 2), 
it is such that 

< d£ >= 0, and < d£ 2 >= 2Vdt. (2) 

Its non-differentiability is apparent here since < 
d£ 2 > 1 / 2 /dt cx dt^ 1 ! 2 and it is therefore formally infi- 
nite. However, as we shall see, the fact that it vanishes 



in the mean allows one to give a description in terms of 
only the differentiable part dx of the elementary dis- 
placement, that includes the indirect effects involved 
by the existence of the non-differentiable part. The co- 
efficient I? is a measure of the amplitude of the fractal 
fluctuations. 

(ii) The test-particles can follow an infinity of pos- 
sible trajectories: this leads one to jump to a non- 
deterministic, fluid-like description, in terms of the 
"classical part" of the velocity field of the family of 
trajectories, v = v(x(t),t). 

(iii) The reflection invariance under the transforma- 
tion (dt <-> —dt) is broken as a consequence of the non- 
differentiability Indeed, there are two definitions of 
velocity which are based on the variation of the posi- 
tion variable X either in the interval [t — dt,t], or in 
the interval [t,t + dt]. The two definitions are equiv- 
alent in the differentiable case, but no longer in the 
non-differentiable one considered here. This leads to a 
two- valuedness of the velocity vector, (u_,t> + ). The use 
of a complex velocity, V = (v + + V-)/2 — i(v + — i>_)/2 
to deal with this two-valuedness can be shown to be a 
covariant and simplifying representation ( Nottale |2()j , 
Celerier 0). 

These three effects can be combined to construct a 
complex time derivative operator ((20]) |2I| . |22|L that 
reads, 

%- = 4- + V.V - WA (3) 
dt dt w 

Now, we can use this new relation and follow the 
standard classical mechanics construction. A general 
Lagrange function L(x, V,t) characterizes any system 
and we define a general action S thanks to 



L(x,V,t)dt. 



(4) 



Since the velocity is now complex, the same is true of 
the Lagrange function and therefore of the action. The 
generalization to a complex velocity V does not modify 
the form of the classical Euler-Lagrange equation (|21j. 
|22|L but now, the classical differential operator d/dt 
is replaced by the complex differential operator d/dt 



d dL dL 
dt dVi dx, 



= 0. 



(5) 



Let /i be the mass of the particle. In the case of 
a Newtonian closed system in a scalar field which is 
considered here, the Lagrange function becomes L = 
■i/iV 2 — This leads to a motion equation that con- 
serves the form of Newton's fundamental equation of 
dynamics 

d. 



»dt V 



-V$, 



(6) 
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which is now written with complex time derivative op- 
erator and complex velocity. The complex velocity V 
still depends on the complex action according to the 
standard form (/xV = V«S). Let us now make a change 
of variable, and define the function: 



1 = 6XP (2^) 



(7) 



In terms of this function, the fundamental equation of 
dynamics now writes 



2i[iV^-(y In /) = V$. 



(8) 



By replacing in this equation the complex derivative 
operator by its expression (Eq. 03), one finds that it 
becomes 



2PV 



4 In/ 

dt 



V 



A/ 
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A 6 



(9) 



The real and the imaginary part of this equation can 
be separated. Indeed, the complex action is rewritten 
S = S + iS , and so, the function (/) takes another 
form: 



/ = exp 



with S the classical action which is connected to the 
classical velocity: fj,V — VS. Finally, the real and 
imaginary part of the motion equation induce the fol- 
lowing system of equations 



p(jL + V . V ) V = _ V (0 + Q), 



dp 

at 



div(pV) = 0. 



(11) 



The first one is a standard Euler-Newton equation but 
with the appearance of an additional potential energy 
term Q that writes: 



Q = -2\iD 



(13) 



and which is a manifestation of the fractal geome- 
try. The second one is the continuity equation. We 
have therefore recovered the current minimal equations 
which are used in the standard hydrodynamical de- 
scription, except for the new potential term. In the 
limit where the fractal fluctuations vanish, i.e. T> — ► 0, 
the potential energy Q vanishes and the description is 
exactly reduced to the standard one. But conversely, 
the introduction of this new potential energy allows 
one to integrate the hydrodynamics equation under the 



form of Eq. [§1 which has the advantage to be linear 
and to allow the obtaining of analytical solutions. 

The generalized hydrodynamical system obtained 
above can now be used as motion equation for a large 
class of systems |S] , namely, all those coming under the 
three conditions that underlie its demonstration: large 
number of possible trajectories, fractal dimension 2 of 
trajectories, and local irreversibility. Actually these 
conditions amount to a loss of information about an- 
gles, position and time. But, paradoxically, such a total 
loss of information on the individual trajectories results 
in a tendency for self organization and structuration of 
the systems [22] • 

3 Theoretical development: 
spherical potentials 

3.1 Symmetry constraints 

The resolution of physical problems depends on the 
conservation laws imposed by the symmetries to which 
the physical system under consideration is subjected. 
As a first step, we shall consider spherically symmet- 
ric systems. Then, for a more complete study, other 
symmetries will be explored, in particular axial sym- 
metrical systems. Moreover, the elementary approxi- 
mation about the spherical potentials will be enriched 
by taking into account dynamical perturbations. These 
various effects are accounted for in a simple way by a 
specific choice of the coordinate systems (e.g. elliptic, 
parabolic or cylindrical). 

Let us consider a test particle ejected by a star, sub- 
mitted to various forces described by potentials (y>i). 
Its motion is described by the hydrodynamical equa- 
tion EaUTl 

M (A+y.V)V = -V(5>i + Q), (14) 

Such a system is expected to be subjected to gravita- 
tional, radiative pressure, electromagnetic, and colli- 
sional forces. 

However, in order to simplify the description we shall 
base ourselves on a fundamental observational result 
concerning these objects: the observations of the shells 
of many PNe's have shown that the expansion veloc- 
ity is nearly constant. Such a result is already used 
in many numerical simulations (|J], ^H], |15|L This 
means that, in the first order approximation (neglect- 
ing the interaction with the ISM), there is a dynamical 
compensation between the various forces acting on the 
shell, in particular between the gravitational force and 
the radiative stellar pressure. We shall therefore as- 
sume, as a first simplifying step, that the potential is 
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almost constant. We take the constant velocity ap- 
proximation (i.e. free test particles) as a starting point 
for the following developments: In a second step, the 
model will include elements of perturbation such as ve- 
locity expansion to second order and the first order of 
perturbation power-series expansion (magnetics contri- 
butions and interactions with the ISM). 

Therefore, we consider the generalized hydrodynam- 
ical equation with a vanishing potential: 



p(^ + V-V)V = -V(Q), 



(15) 



The velocity seems to be invariant during the ejection 
process. To simplify the search of stationary solutions, 
it is easier to rewrite the Euler-Newton/continuity sys- 
tem (Eq. 1111 Eq. ^] depending on the variables [p, V]) 
in terms of a unique complex equation (depending on 
the variables [r, t]) since it takes a Schrodinger-like 
form, for which several analytical solutions are known. 
This equation emerges from the integration of Eq|5]and 
writes f|20]. 



DA/ 



df $ 



0. 



(16) 



In a constant potential, stationary solutions of the form 
f(r, t) — g(r) exp(—iEt/2pD), can be searched, and we 
obtain: 

2pV 2 A g(r) - Eg(r) = 0, (17) 

where E = p 2 /2p — 2pD 2 k 2 is the energy of this free 
particle. If we take a constant potential different from 
zero, we obtain the same equation but with a differ- 
ent particle momentum. The property of the Lapla- 
cian in the spherical coordinates governs the use of the 
spherical harmonics (which arc common functions of 
L 2 and L z ). The general solution is now constrained 
by the conservation of the angular momentum and we 
can take a solution of the form g(r) = R{r) Y[ a (9 1 <fi), 
where I and m are integers. The radial part is solution 
of the equation 



R (r) 



-R (r) + 
r 



,2 1(1 + 1) 



k Z - 



R(r) = 0. (18) 



Two radial solutions of this free particle equation are 
known |17|: 



1 d \ 1 e ±ikr 



(19) 



These functions can be developed in terms of the first 
order Hankcl functions, namely 

fl±(r) = ±iA\I^H%l\kr). (20) 



which represent divergent spherical waves (+ case, with 
first order Hankel function) or convergent spherical 
waves (— case, with second order Hankel function). 
Divergent spherical waves can correspond to a flow of 
central particles emission and the convergent one to an 
infall flow of particles (we can already note that the 
physical solutions are totally reversible: this particular 
behavior will be developed in a forthcoming section). 
The normalization condition restricts the value of the 
constant, so that A = l/V2Vk. The square of the 
modulus of the function / gives the probability den- 
sity: 

| 5 (r)| 2 = |Jfa(r)| 2 |ir(M)| 2 . (21) 
We expect the true matter distribution to follow pref- 
erentially the peaks of this probability density distri- 
bution. 

Now, this equation represents the distribution of po- 
tential trajectories for particles ejected during a unit 
time. But our aim is to know the evolution of these 
trajectories for distances and times higher than the 
ejection area. As a first step, the emitted particle flow 
can be considered to be only radial. Thus, the pro- 
gression of the flow is free along the radial coordinate, 
and only the radial part of the function depends on the 
time variable. Thus the progression of the maximum 
probability can be approximately described as a spher- 
ical wave front with a constant velocity Vq. The time 
dependent function simply reads: 

g(r,t) = ^R kl {r-V t)Yr{6,4>). (22) 

The 1/r 2 coefficient accounts for the dilatation of the 
elementary shell. 

3.2 Effect of perturbations 

Let us now go on with our progressive method and 
include, for a more complete description, the pertur- 
bative potential terms. 

First, one can relax the hypothesis of constant ve- 
locity by the consideration of a power series expansion 
of the radial velocity. These radial perturbative terms 
only affect the radial part of the spherical potential. 
Then this restriction will only change the radial part 
of the solution, but not the angular part. The infor- 
mation about the matter density distribution along the 
shells is related to the ejection history, while the global 
shape (which is observed in the end) is given by the 
angular distribution of matter (theoretically described 
in our framework by spherical harmonics). Thus, the 
principal actors of the structuration are the symmetries 
of the system. We can generalize this argument by in- 
cluding all kind of radial perturbative potentials (com- 
patible with ejection processes) without any change in 
the result obtained about morphogenesis. 



4.2 Comparison with observational data 
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Other perturbations could affect the global process. 
We describe, hereafter, the effect of two of them: 

• Self gravitation: The test-particle ejection associated 
with bipolar solutions (initially associated with bipolar 
axisymmetrical spherical harmonics) evolves through a 
new dynamical perturbation field linked to the axial 
symmetry. A self gravitational potential is created by 
the initial cylindrical repartition: each particle is sub- 
jected to a perturbative force in — 1/p 2 (where p is the 
radial cylindrical coordinate). This attractive gravi- 
tational force constricts the ejection cone toward the 
symmetry axis of the system. 

• Magnetic field: A magnetic contribution must be in- 
troduced in terms of a perturbative field. Strong mag- 
netic fields are naturally associated with several stellar 
stages (e.g. red-giant stars and white dwarf stars [T%|)- 
A poloidal field can be introduced, as a first approxi- 
mation, and the evolution of ionized particles (because 
of the UV radiation) is described by the same dynam- 
ics: a constriction along the axial system symmetry. 

4 Results and comparisons to 
the observational data 

4.1 Theoretical results 

The main result obtained in the present work is the 
expected discretization of the possible shapes and the 
quantization of conservative quantities associated to 
the various predicted morphologies. 

Let us first consider the main conservative quantity, 
i.e., energy. The energy to mass ratio (E/p = 2V 2 k 2 ) 
is not quantized in this case, since there is no bound- 
ary condition to restrict the value of the particle wave 
vector (k). However, one will be led in future devel- 
opments to consider the connection between the initial 
ejection process and the internal star structure. This 
work could provide a restriction of the wave vector pos- 
sible values. 

The second and the third prime integrals (consis- 
tent with the spherical coordinates) involve a quanti- 
zation of the square of the reduced angular momentum 
and of its projection on the z axis, namely, (L/p) 2 = 
21(1 + 1)V 2 and L z /p = 2mT>, which are independent 
of the test mass. The ratio L/L z = y/l(l + l)/m is also 
independent of the scale-parameter D. Thus, this rela- 
tion gives a general description of the angular momen- 
tum conservation law in all spherical ejection processes, 
independently of the specificity of the underlying frac- 
tal process. 

The interpretation of these solutions is that the mat- 
ter is expected to preferentially fill the high probabil- 
ity regions, as defined by the geometry of the geodesies 



(provided there is a sufficient quantity of ejected mat- 
ter). Thus, one can directly associate a quantized 
structure with the angular functions. Moreover, as we 
shall now see, spherical harmonics involve also a quan- 
tization of other observables (i.e., the existence of peaks 
in their probability density distribution). 

The 9 angle is also subjected to a probability den- 
sity distribution that have quantized peaks. This an- 
gle is observable in terms of the initial ejection an- 
gle (in the nearest central area, where the pertur- 
bations are insignificant). Concerning the remaining 
structure, the theoretically predicted value of the an- 
gle is no longer observable since the perturbation ef- 
fects (studied above) result in a deformation of the 
shape toward the symmetry axis of the system. An- 
other consequence of using spherical harmonics for the 
morphology description is the dependence of the mat- 
ter density distribution on the angle quantization. For 
each component of the system (i.e. axial ejections, 
bipolar shells or equatorial discs), one can predict an 
initial ejection angle probability and a relative matter 
density distribution (relative to the other components 
of the matter distribution). All these results about 
the expected more probable values of the observable 
variables are summarized in Tabled F°r each couple 
(l,m), we give the most probable values of the ratios 
L 2 /p, L z /p and we plot the associated angular solu- 
tions, p(9) — |YJ m (#, csi)| (in a simplified form that 
takes account of the symmetry under the <fi rotation) . 

Another important information in this table con- 
cerns the ratios of the amplitudes of the probability 
peaks for the ejection angles. The relative matter den- 
sity between the imbricated conic structures (see Table 
EJ is expected to be proportional to the relative prob- 
ability density: such a theoretical expectation can be 
put to the test (by an analysis of the luminosity ratios) 
in the observational data (see example hereafter). 

The various shapes allowed by the generalized hy- 
drodynamical solutions are presented in Table [5] The 
bipolar shells arc plotted with a weak deviation to- 
ward the axis of symmetry. This deviation is typical of 
the various perturbative effects considered here-above 
(e.g., self gravitation, magnetic field). 

4.2 Comparison with observational 
data 

4.2.1 Planetary nebulae 

As recalled above, the allowed shapes only depend on 
the (I, m) couple. Thus, the expected morphologies 
are no longer a consequence of the only initial condi- 
tions; they are mainly linked to the values of the con- 
servative quantities and related to the angular bound- 
ary conditions. Similar basic hypotheses allow one 
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Table 1: Synthesis of the quantized values of the squared angular momentum L 2 /fi, of its projection along z axis L z /fj,, and of 
the matter distribution along quantized angles p(8), associated with the angular wave functions Y™, for the I and m = 1 to 7. 
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Table 2: Synthetic representation of the first quantized shapes resulting from the angular ejections solutions given in (Table HI . 
Note the analogy of structure of the shapes along diagonals, i.e. I — m = k =cst, which are characterized by the same number of 
imbricated cones, and differ only by the values of the angles. Note also that the theoretically predicted density significantly varies 
between the cones (see the corresponding figures in Table fT\ . 
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Figure 1: The elementary isotropic structure emerges 
from the fundamental level of the angular quantization 
law (I = 0, to = 0). Left hand side: theoretical predic- 
tion. Right hand side: the spherical planetary nebula 
IC 3568. 



Figure 2: For some values of the inclination, the theo- 
retical equatorial disc (left hand side) could appear like 
an elliptic PN. For example, NGC 3132 is observed el- 
liptic, but there is still discussions about its true mor- 
phology (Monteiro et al. 2000). 



to understand the constitution of bright shells in the 
spherical, elliptic and bipolar planetary nebulae. For 
remaining uncommon morphologies, various perturba- 
tions could be included in the active potentials. More- 
over, more complicated symmetries (linked to other dy- 
namical conservation laws) remain to be explored. 
A new morphological description of planetary 
nebulae: 

• Spherical: For a (l,m) couple equal to (0,0), one ex- 
pects an isotropic matter distribution. This elementary 
case reflects the properties of the spherical planetary 
nebulae which were first discovered. 

• Elliptic: All the solutions consistent with (I = m) 
involve an equatorial disc structuration. With differ- 
ent tilting on the line of sight, elliptic shapes emerge 
from this solution. It is important to note that el- 
liptic shapes can also result from perturbed spherical 
shapes, in accordance with the ISW model for elliptic 
shapes. Moreover, provided some apparently elliptic 
shapes are indeed the result of such a projection effect 
of disk structures, one expects to observe one day ex- 
tremely flattened elliptical PNe: however, up to now 
this kind of PNe have not yet been observed (to our 
knowledge). 

• Bipolar: Bipolar shapes emerge naturally from this 
work. These shapes are, in the standard approach, ex- 
plained in terms of various effects such as different den- 
sity outflows, presence of a double star system, etc... 
The present approach does not contradict these views: 
we simply connect the shapes directly to the state of 
the system, as described in terms of its conservative 
quantities (at first, the state of angular momentum). 
Whatever the particular cause for the system to jump 
to this state (which can be the presence of a second 
star, etc.), the same shape is expected to be observed 
for a given state. Moreover, this result generalizes the 
ISW model principles that lead to bipolar structures, 




f 

HbT2 



Figure 3: Elementary bipolar structure: the simple 
shell structure. The predicted (5,4) quantization with 
an initial angle of 63° is consistent with the observed 
structure of the planetary nebula Hbl2 . 



since the free particles in bright shells could be directly 
issued from wind shocks (which are not conserved in 
standard models). 

The first bipolar structure allowed by the possible 
values of the angular momentum is the simple shell 
one. This kind of simple situation can be illustrated 
by the planetary nebula Hbl2 (fig. 0J. The initial an- 
gle is consistent with a (I = 5, to = 4) quantization 
(near 63°) but this value is not sufficiently well mea- 
sured to specify the quantum numbers with certainty 
(with regard to the nearest possibilities: (I = 4, to = 3) 
or (I = 6, to = 5)). 

Bipolar shapes become more complicated when the 
angular momentum increases. Beyond simple shells 
(to = I — 1), bipolar shells with equatorial disks emerge 
in the same way (to = 1 — 2, cf thciJ65 structure). 

The next case (fig. 0} is the double-shell plane- 
tary nebula of the M2-9 type. This shape is naturally 
expected from our classification, as corresponding to 
to = I - 3 (table EJ). The morphology of M2-9 can 
be associated to the quantum numbers (5, 2) (cf. ta- 



4.2 Comparison with observational data 
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Figure 4: Double shell structure: the M2-9 planetary 
nebula can be associated with an angular momentum 
quantization of quantum numbers (I = 5, m = 2). 

ble^l which involves a quantized angular distribution 
around the values (33°, 72°, 108°, 147°) and a higher 
density for the inner shells. The PN structure near 
the central star is not clearly defined, but the extrap- 
olation of the external shell allows us to estimate the 
values of the initial ejection angles. This estimation 
(ss 30° for the outer shell and « 70° for inner one) 
is consistent with the predicted angles: moreover, the 
inner observed density is clearly higher than the outer 
shell density, as theoretically expected. 
Even more complex (therefore less probable) structures 
that have not been observed up to now, are still pos- 
sible: double-shell + disk (to = I — 4), triple-shell 
(m = I — 5), etc... In all these new cases, as for the al- 
ready observed ones, the theoretical prediction applies 
not only to the overall morphology, but also to quan- 
titative values of angles, angular momenta and density 
distributions, and it will therefore be possible to put it 
to the test in the near future. 

4.2.2 LBV stars 

The theoretical development presented in the first part 
of this paper is associated with a general ejection 
process. Luminous blue variable (LBV) stars are in 
essence very different from PN central stars. Despite 
of this internal incompatibility, some recent works |TT] 
have been developed to explain the shell formation 
around LBV stars by the same hydrodynamical ISW 
model that was used for planetary nebulae. 

Therefore the method studied in the present paper 
can be applied to a morphological description of these 
objects. In consequence we also expect a quantization 
of quantities such as angular momenta and angles and 
a discretization of the possible shapes for the matter 
density distribution around LBV stars. For example, 
the Eta Carinae LBV star shows bipolar shell with an 
equatorial disk (fig. such a structure is expected in 
all the cases when I > 3 and to ~ (1 — 2)). There- 
fore this result supports the idea of a general behavior 



Figure 5: Generalization to another singular stellar sys- 
tems: comparison between theoretical shapes with sim- 
ple bipolar shell/equatorial disc and a LBV star (Eta 
Carinae) . 

of ejection processes and of the self-similarity of shell 
structurations, independently of the nature of the cen- 
tral star. 

4.2.3 Supernovae 

There are other systems that can be analyzed with 
this new point of view, in particular supernova rem- 
nants. Whatever the specific ejection process at work 
in supernovae (la or II..), the three conditions upon 
which our method relies, namely (i) fractal trajecto- 
ries, (ii) infinity of potential trajectories, (iii) time- 
reflection breaking, are expected to be achieved in this 
case, since the same collisional/diffusing processes as 
for planetary nebulae are present. Moreover, the here- 
above theoretical developments can be generalized to 
scattering systems described by an initial convergent 
wave function (as in implosion processes) and by a re- 
sultant divergent wave function (ejection processes). 
Supernova structures emerge from the interaction of 
remnant matter from the central star explosion and 
the ISM. This simple presentation suggests a behavior 
of SNe particles which may be similar to ejected par- 
ticles in PNe or LBV stars. Thus, provided the ISM 
interaction is assumed to be homogenous (as a first 
order approximation), the quantized results obtained 
by the theoretical development (table |2J are still valid. 
The account of inhomogeneities of the ISM would be 
more difficult, since they degrade the initial conserved 
information on prime integrals and then, would involve 
a more complicated description. One of the most enig- 
matic SN is SN-1987A (Fig. EJ). Its particular structure 
(which has up to now resisted to all hydrodynamical 
simulations) can be accounted for by a unique ejection 
associated with a quantized angular momentum corre- 
sponding to to = I — 2 (see Table|2Jl. This configuration 
induces, in the case of a discontinuous ejection process, 
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Figure 6: Simple model of supernova remnant: as a 
first order approximation, one finds that supernovae 
remnants could also be structured like other spherical 
ejection systems. For a single ejection phase, this leads 
to the formation of expanding rings, that are expected 
to be ejected with higher probability along given an- 
gles. 

a geodesic structuration into two bipolar rings and an 
equatorial ring, as observed. Moreover, in recent im- 
ages the central star is projected on the right ring. In 
such a configuration, the ejection cone angle is given by 
cos(#) — b/a where a and b are respectively the major 
and minor axes of the projected ellipse. From 6 mea- 
surements on the inner edge, middle and outer edge of 
the equatorial ring and of the right polar ring (the left 
ring was not used because it shows distortions to ellip- 
ticity), we find 9 = 41.2 ± 1.0 deg, which supports the 
identification with the case m = 2, 1 = 4, for which the 
predicted angle is 40.89 deg (while the nearest similar 
configurations yield 31 deg and 47 deg). The youth of 
the system (10 years compared with 10 4 years for typ- 
ical PNe) supports the assumption that the observed 
angle is still the initial ejection angle. 

From this preliminary work, we suggest that future 
observations of young supernovae in homogeneous ISM 
will allow one to put the theory to the test in a quanti- 
tative way, since we expect the occurrence of quantized 
structures consistent with the general results given in 
Table 

4.2.4 Star formation 

Low mass star formation processes seem to be univer- 
sally associated with circumstellar infall / accretion 
discs and outflow / axial ejections (pQ, [S]). In the 
standard approach, the existence of such simultaneous 
structures is still not fully understood (despite several 
attempts using various magnetic models [T2|1. 

The solutions associated with free motion in spheri- 
cal potentials describes either convergent or divergent 
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Figure 7: An example of morphology observed in a 
star formation process, compared with a theoretically 
predicted possible structure. 

spherical waves (eq. I19JI . From the point of view of 
the present theoretical description in terms of probabil- 
ity amplitudes, the two cases of ejection and accretion 
both correspond to a scattering process (one simply re- 
verses the sign of time). It is therefore quite simple to 
consider solutions that combine accretion (for example, 
equatorial) and ejection (fro example, polar). 

One cannot make any longer the hypothesis of free 
motion in an accretion phase. However, this concerns 
only the radial dependence of the solutions, while the 
angular dependence is still expected to be given by 
spherical harmonics. Finally, in a typical star birth 
stage, the above method is able to account for the ob- 
served discretized structuration of the infall matter and 
of the resultant particle ejection. 

For example the HH-30 observations (Fig. [JJ allow 
us to make a preliminary analysis of this proposal. The 
equatorial disc and the resultant axial ejection agree 
with the morphology given by I = 1, m = for axial 
ejection and m = I for the accretion disc. The global 
angular momentum is conserved during the scattering 
process, but the L z behavior is still an open problem 
(since it could be described by a more complex solu- 
tion like (I = 2,m = 0)). The bright bipolar shells 
are currently considered as reflection nebulae pi], but 
this component could also have a different dynamical 
behavior (with respect to the equatorial disc) and thus 
could correspond to a solution like (3,0) or (4,0). We 
present in Fig. (0 the structure closest to observations 
allowed by the theoretical developments. 

5 Conclusion and prospect 

In this paper, we have given a general description of 
matter structuration in spherically symmetric poten- 
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tials. The motion equation, rewritten in terms of an 
hydrodynamical system which incorporates the equa- 
tion of continuity, depends directly on the conditions 
of symmetry of the system. 

In the case of spherical symmetry, which is the sit- 
uation mainly considered in this study, the solutions 
depend on the values of the conservative quantities (E, 
I? and L z ). This method leads to a great diversity of 
stationary solutions, that are expected to describe in- 
fall and ejection motion from central bodies, in terms of 
probability density distributions for the radial distance 
and the values of the ejection (or accretion) angles. 

In the case when the potential depends only on the 
radial coordinate, the general solutions (which give the 
density of potential trajectories) are then separable in 
a radial R(r) and a spherical part given by the spheri- 
cal harmonics Y l m (8, 4>). The possible morphologies are 
therefore classified according to the values of the inte- 
ger numbers I and m. Thus, most probable quantized 
structures (table EJ result from the angular probabil- 
ity distributions (table and can be compared with 
the various systems that come under such a descrip- 
tion (planetary nebulae, LBV stars, supernovae, young 
star formation). Moreover, the theoretical approach al- 
lows one to predict the possible existence of structures 
that have no yet been observed, and of yet unidentified 
correlations between the morphology and observable 
quantities. 

Let us conclude by some examples of possible future 
extensions of the present study: 

- The connection between internal stellar structures 
and the ejection process could provide informations 
about the probability distribution of the velocity ejec- 
tion. 

- Concerning the planetary nebula description, the 
internal stellar angular momentum is known to be 
structured: various external quantized solutions are 
then possible. The internal structuration could induce 
an external combination of spherical harmonic solu- 
tions (depending on the ejection history ). 

- Particular quantized solutions (l,m) can be con- 
nected with the central system and show a dependence 
on specific conditions (single star with inclined mag- 
netic fields, double star system etc.). 

- Several simple systems present axial symmetry: 
therefore several new possibilities of morphogenesis can 
be explored by considering other coordinate systems 
(e.g. cylindrical symmetry, parabolic symmetry etc.). 

- Only hints about the perturbation effects have been 
given here: future work is needed in order to introduce 
the perturbation potential in the Euler-Newton / con- 
tinuity equation (non trivial solutions). 

- Numerical simulations using the above solutions 
as initial conditions can be performed, in order to 



study the influence of perturbation fields (e.g. mag- 
netic contribution, self-gravitation, ISM interaction...), 
and to explore the possibility of appearance of non- 
trivial structures. 
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